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Detection of anhydrous hydrochloric acid, HCI, in IRC+10216 with 
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ABSTRACT 

We report on the detection of anhydrous hydrochloric acid (hydrogen chlorine, HCI) in the carbon-rich star IRC+10216 using the spectroscopic 
facilities onboard the Herschel satellite. Lines from J= 1-0 up to J=7-6 have been detected. From the observed intensities, we conclude that HCI is 
produced in the innermost layers of the circumstellar envelope with an abundance relative to H 2 of 5xl0~ 8 and extends until the molecules reach 
its photodissociation zone. Upper limits to the column densities of A1H, MgH, CaH, CuH, KH, NaH, FeH, and other diatomic hydrides have also 
been obtained. 
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1. Introduction 

The chemistry of chlorine (CI) in the interstellar medium (ISM) 
is particularly poorly known, mainly because of the relatively 
small number of C l-containing molecules detected to date (see 
Neufeld & Wolfire 2009). This element has two stable isotopes 
( 35 C1 and 37 C1) and h as a relatively l ow so lar abundance of 
3xl0" 7 , relative to H dAsplund et al. L |2009|) . Anhydrous hy- 
drochloric acid (HCI; also known as hydrogen chloride) remains 
the only chlorine-bearing molecule observed to date in the inter- 
stellar medium, and is believed to be one of the major reservoirs 
of chlorine in the ISM. It ha s been observ ed in both th e dense and 
diffuse interstellar medium dBlake et al . , 1985: Feder man et al. L 
1 1995b . The three metal chlorides A1C1, NaCl, and KC 1 have also 
been observed in space (Cernicharo & Guelin , 1987), but solely 
in circumstellar envelopes (CSEs) around evolved stars, where 
they are formed in the hot and dense stellar atmospheres under 
thermochemical equilibrium. Although it has not yet been ob- 
served in such environments HCI is also expected to be a ma- 
jor chlorine species in circumstellar envelopes. Being a light hy- 
dride, its rotational transitions lie in the submillimeter and far- 
infrared domain, which is difficult to observe from the ground 
because of severe atmospheric absorption. Other light species, 
mostly metal-bearing hydrides such as A1H, FeH, Mg H, and 
CaH, are detected in sunspots and M-type stars dGizis L Il997t 
Wall ace et al. L 1200 ll) and are also potentially present in the in- 
nermost zones of carbon-rich circumstellar envelopes. 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA 



The infrared source IRC+10216 (CW Leo) is one of the 
brightest in the sky, making it an ideal target to be observed with 
the Herschel Space Observatory (Pilbratt et al. , 2010). Around 
50% of the molecules observed in space have been detected to- 
wards this object. Most of those molecules are heavy carbon 
chain r adicals dCe rnicharo & Guelin , 1996a), metal-bearing 
species dCernic haro & Guelin , 1987), and bo th diatomic and 
triatomic molecules (Cernicharo et al. L l2000h . Its far-infrared 
spectrum, obtained with low spectral resolution with the Infrared 
Space Observatory (ISO), was analyzed by Cernic haro et al. I 
(Il996bl) with the Infrared Space Observatory (ISO) with lim- 
ited spectral resoluti on. The spectrometers on board Herschel 
(Pilb ratt et al. L l2010l) offer the possibility to search for light di- 
atomic hydrides with high sensitivity, thanks to the telescope's 
large collecting area and the performances of the instruments, 
and with high spectral resolution compared to ISO. In this Letter, 
we report on the first detection of HCI toward the circumstellar 
envelope of the carbon-rich star IRC +10216, and discuss the 
implications for the chemistry of chlorine in these astronomical 
regions. We also present upper limits to the abundance of metal 
hydrides. 



2. Observations and data reduction 

The three instrume nts on board the Herschel satellite 
(Pilb ratt et al.L 1201 Oh have medium to high spectral resolu- 
tion spectrometers. PACS and SPIRE spectroscopic observa- 
tions were obtained in the context of the guaranteed time 
Key programme Mass-loss of Evolved StarS (Groenewegen et 
al., in prep). The PACS instrument, its in-orbit performance 
and calibration, and its scientific capabilities are described in 



J. Cernicharo et al.: Detection of HC1 in IRC+10216 



Pogli tsch et al.. I (120101) . The PACS spectroscopic observations 
of IRC +10216 consist of full SED scans between 52 and 210 
/im obtained in a 3x1 raster, i.e., a pointing on the central ob- 
ject, and two additional ones located 30" on each side of it. The 
observations were performed on Nov 12 2009 (OD 182). The po- 
sition angle was 110 degrees. The instrument mode was a non- 
standard version of the chop-nod PACS-SED AOT, used with a 
large chopper throw (6'). A descriptio n of that mode and o f the 
data reduction process can be found in Rove r et al.. I ([2010). The 
estimated global uncertainty in the line fluxes is 50 %. However, 
the relative calibration is much better, hence it is possible to esti- 
mate the contribution of the most abundant species to the lines of 
HO by using adjacent lines of these species. PACS and SPIRE 
photometry observations are presented in Ladjal et al. (2010). 

The SPIRE FTS measures the Fourier transform of the 
source spectrum across short (SSW, 194-313 //m) and long 
(SLW, 303-67 lyum) wavelength bands simultaneously. The 
FWHM beamwidths of the central SSW and SLW pixels vary 
between 17-19" and 29-42", respectively. The source spectrum, 
including the continuum, is restored by taking the inverse trans- 
form of the observed interferogram. For m ore details about the 
SPIRE FT S and its calibrat ion, we refer to iGriffin et all d201dh 
and Swin vard et al.. I d2010h . We made use of two observations 
of IRC+ 1 02 1 6 with the high-resolution mode of the SPIRE FTS 
on the 19 November 2009 (OD 189). For each observation, ten 
repetitions were carried out, each of which consisted of one for- 
ward and one reverse scan of the FTS, each scan taking 66.6s. 
The total on-source integration time for each FTS spectrum of 
IRC +10216 was 1332 seconds. In the end, both FTS spectra 
were averaged. The unapodized spectral resolution is 1.2 GHz 
(0.04 cm"'), which after apodization (using extended Norton- 
Beer function 1.5; Naylor & Tahic 2007) became 2.1 GHz (0.07 
cm -1 ). The sensitivity of the SPIRE/FTS spectrometer allows us 
to detect lines as weak as 1-2 Jy. The whole PACS+SPIRE spec- 
trum of IRC+10216 has been shown by Decin et al. (2010). 

3. Results 

In addition to lines arising from v ibrational levels up to 10000K 
dCernicharo et al. L Il996bl l2010h . HCN is the main contribu- 
tor to the spectral features detected with the SPIRE and PACS 
spectrometers. Th e frequencies of these lines were calculated by 
Cernicharo et al. (2010) from the rotational constants provided 
by iMaki et aLl (1199 6. 2000). These frequencies were used to 
identify most features shown in Fig. 1. The other strong features 
arise from CO, SiS , SiO, and CS. These lines were analyzed by 
iDecin et al. I (J2010), and the data used to estimate the contribu- 
tion of these species to the lines of HC1 shown in Fig. 1. HC1 has 
two stable isotopologs, H 35 C1 and H 37 C1, the former being 3.1 
times more abundant than the latter in IRC +102016, accord- 
ing to the 35 C1/ 37 C1 abundance ratio derived from observations 
of NaCl, KC1, and A1C1 dCernicharo et al. L |2000T) . Frequencies 
for H 35 C1 and H 3 7 C1 were computed f r om th e rotational con- 
stants derived by ICazzoli & Puzzarini (2004). The laboratory 
measurements have an accuracy better than 0.5 MHz for lines 
up to J,,^,— 14. 

The first two rotational transitions of HC1 are covered within 
the spectral range of SPIRE, and are detected as emission lines in 
the spectral data obtained towards IRC +10216 (see Fig.Q}. The 
low spectral resolution of SPIRE (2.1 GHz) prevents us from 
resolving the individual emission components related to H 35 C1 
and H 37 C1, which are separated by 940 and 1879 MHz for the 
7=1-0 and 7=2-1 transitions, respectively. The 7=1-0 emission 
is blended with the 7=13-12 transition of C 34 S. Together, they 



appear as a shoulder at the high frequency side of a stronger 
emission feature, which is a composite of several rotational lines 
of HCN in different vibrational states (we estimate a flux of 3 Jy 
for the J=l-0 line of HC1). The 7=2-1 transition is less severely 
blended with stronger features, although it does overlap with the 
7=26-25 transition of C 34 S and with some ^-doubling compo- 
nents of the 7=14-13 transition of HCN in its V2=4 vibrational 
state. Inspection of the line intensities of these species (C 34 S and 
HCN V2=4) in the nearby spectral region (see Decin et al., 2010), 
indicates that they contribute at most half of the intensity of the 
detected emission feature. Hence, the measured flux for the 7=2- 
1 transition of HC1 is 9 Jy. 

The spectral range of PACS covers the J=3-2 up to the J=9- 
8 rotational transitions of HC1. The data acquired toward IRC 
+ 10216 allow one to clearly identify the 7=3-2 to 7=7-6 tran- 
sitions (see Fig. [TJ, the components related to H 35 C1 and H 37 C1 
being spectrally resolved. The lines are clearly seen without the 
contamination of stronger lines with derived fluxes, after remov- 
ing the contribution from other species, of 12, 26, 40, 45, and 
60 Jy for the J=3-2 to J=7-6 lines, respectively. Among them, 
only the 7=4-3 rotational transition appears appreciably con- 
taminated, mostly by the 7=29-28 transition of H 13 CN, which 
severely hampers the visualization of the H 37 C1 line, and to a 
lesser extent by the 7=28-27 transition of HCN in its V2=3 t=\ 
vibrational level, which still leaves the H 35 C1 line visible. The 
H 35 C1 component of the J '=7-6 transition is observed as an emis- 
sion feature at the correct frequency and with an intensity com- 
patible with the lower 7 transitions. However, the signal-to-noise 
ratio is only 5, providing only upper limits to the intensity of the 
H 37 C1 isotopomer. Higher J lines are not detected because of the 
limited sensitivity at these frequencies. In spite of the low spec- 
tral resolution and the, in some cases quite severe, contamination 
of the observed HC1 lines, the large number of transitions cov- 
ered by the SPIRE and PACS data makes the identification of 
HC1 in the circumstellar gas of IRC+10216 quite certain. 

The observed HC1 lines have been interpreted with the aid of 
an excitation and radiative transfer model based on the Large 
Velocity Gradient (LVG) formalism. The H 35 C1/H 37 C1 abun- 
dance ratio is poorly constrained from the observations and 
was fixed to be 3.1, as derived for the 35 C1/ 37 C1 abundance ra- 
tio from previous observations of NaCl, KC1, and A1C1 in IRC 
+ 10216 dCernicharo & Guelin L fl987t ICernicharo et al. L 12000). 
We included the first 20 rotational levels within the ground vi- 
brational state of bot h H 35 C1 and H 37 C1. The adop ted dipole 
moment is 1.109 D dDe Leeuw & Dvmanus L Il97lh . The rate 
coefficients for collisional de-excitation from HC1 levels up 
to J=7, through collisio ns with H2 and He, were taken from 
iNeufeld&Greenl d!994h . and the Infinite Order Sudden (IOS) 
approximation was applied for higher 7 levels. The circumstellar 
envelope is simulated as a spherically distributed gas expanding 
at a constant velocity of 14.5 km s -1 . Th e gas densi t y and tem- 
perature radia l profi les were taken from lAgundezl d2009l) and 
iFonfria et ail d2008l). Th e adopted distance to IRC +10216 is 
120 pc dSchoier & OlofssonU2001l) . 

Molecules in IRC + 1 02 1 6 are either concentrated around the 
central star (e.g., HCN) or distributed in a hollow shell of radius 
10"-20" (e.g. CN). The lack of information about the HC1 line 
profiles prevents one from deciding which of these two distri- 
butions HC1 follows. Several radial distributions for HC1 were 
tested, but the LVG model indicates that the observed relative 
intensities of the HC1 lines can only be reproduced if HC1 is 
concentrated around the central star. We thus adopted an abun- 
dance radial profile in which the abundance of HC1 relative to 
H2 is constant from the stellar photosphere out to the radius 
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where it is photodissociated by the ambient interstellar UV field. 
The adopted photodissociation rate for HC1 is l.lxlO -9 exp(-1.8 
Ay) s _I (Roberge et al. , 1991), where Ay is the visual extinc- 
tion against interstellar light at each radial position in the enve- 
lope. The derived abundance of H 35 C1, relative to H2, is 5xl(T 8 , 
which produces line profiles that are in reasonable agreement 
with the observed ones, as shown in Fig.Q] 

3.1. Chlorine chemistry in IRC +10216 

To verify whether the above conclusion is compatible with 
chemical arguments, we calculated the composition of the gas 
in the surroundings of the stellar atmosphere of IRC +10216 un- 
der therm ochemical equilibrium (TE). The utilized code is de- 
scribed in iTeiero & Cernicharo I d!991l) . We inclu ded 24 chem- 
ical elements with solar abundances (Asplu nd et al. L [2009), 
and assumed a higher carbon abundance so that [C]/[0]=1.5. 
The t hermochemica l data of the included molecules were taken 
from IChase I d!998l) . with updates for various species being 
taken from the recent literature. In particular, th e thermochem- 
ical d ata of TiH were updated according to Bur rows et al. I 
(200 5J). The adopted parameters for IRC+10216 are presented 
in lAgundez & Cernicharo (2006). In Fig. [2] we show the cal- 
culated abundances of several hydrides in the stellar atmosphere 
and inner envelope of IRC +10216. We see that HO reaches 
a maximum abundance of 6xl0~ 7 relative to H2 (i.e., half that 
value if expressed relative to H) in the 2-3 R* region, where 
TE still holds. According to the results shown in Fig. [2] the 
species that lock most of the chlorine in the inner envelope of 
IRC +10216 are atomic CI (in the region inner to 2 R«), HC1 
(in the 2-3 R» region), and A1C1 (at radii larger than 3 R„). 
The HC1 abundance derived from the SPIRE and PACS data is 
lower than the calculated thermochemical value by a factor of 
~10, which is acceptable given the uncertainties associated with 
the observations and with both the LVG and the thermochem- 
ical equilibrium models. Other Cl-beari ng molecules observed 
in IR C +10216 are A1C1, NaCl, and KC1 (ICernicharo & Guelin L 
1987), for which the abundances relative to H2 derived from the 
IRAM 30- m data are 3.5X 10' 8 , 1.0xl0~ 9 , and 2.5xl0 -10 , re- 
spectively ( Agundez , 2009). Hence, HC1 and A1C1 are the most 
abundant chlorine-bearing molecules in the inner envelope of 
IRC +10216, HC1 being slightly more abundant. As the circum- 
stellar gas expands, molecules become exposed to the ambient 
UV field and are then photodissociated. The chlorine carried 
by HC1 will then be liberated to the gas phase as atomic CI, 
which in the outer molecular shells (at 10" -20") is just slowly 
ionized by interstellar UV photons (the ionization potential of 
CI is 13.0 eV, i.e., lower than, but very close to, that of hydro- 
gen). Therefore, it can participate in rapid chemical reactions to 
form new Cl-bearing molecules. Which are the most likely such 
molecules is difficult to predict because of the lack of chemi- 
cal kinetics data for reactions between atomic CI and abundant 
molecules in the outer envelope of IRC +10216. The literature 
on the chemical kinetics of chlorine is vast, but focuses mainly 
on species such as chlorofluorocarbons (CFCs), which are of in- 
terest as important pollutants of the terrestrial atmosphere. Some 
reactions of atomic CI with unsaturated hydrocarbons such as 
C2H2 have been also studied, although the derived rate constant 
is normally for the three body process, which is not of inter- 
est in the outer env elope of IRC +10216 The UMIST Database 
for Astrochemistry dWoodall et al. LI2007I) contains several reac- 
tions involving chlorine species, but is mostly oriented toward 
the formation of HC1 in interstellar clouds. Since atomic CI is 
relatively reactive, it is expected to react with abundant carbon- 



bearing molecules, such as the radicals C2H or C4H, forming 
chlorine-carbon molecules. The detection of these species (e.g., 
CC1, HCC1, C2CI, etc.) would certainly aid in understanding the 
chemistry of chlorine in the outer layers of circumstellar en- 
velopes, which until to now has been a mystery. 

3.2. Other hydrides in IRC +10216 

As stated above, several metal hydrides have their rotational 
transitions in the submillimeter and far-infrared domains. These 
species are known to be abundant in the photosphere of AGB 
stars dGizislfl997h . and we expect to detect them if they are not 
condensed onto dust grains. We checked the rotational transi- 
tions of several of these hydrides and we found no clear detection 
at the sensitivity limit of the data (3 Jy for SPIRE and 10-20 Jy 
for PACS, depending on the degree of blending, at the 3<x level). 
Moreover, many of their transitions are located in frequency 
domains in which the line confusion limit is reached. This is 
mainly the case for the ranges associated with the HCN lines. 
As an example, the J=5-4 and 10-9 lines of A1H lie at 62.825 
and 124.562 cm -1 , respectively. Neither line is detected (see Fig. 
1). Assuming that the emitting region for metal hydrides has an 
average kinetic temperature of 300 K, then the 3<x upper lim- 
its to their abundance (relative to H 2 ) are 3xl0~ 9 (NaH), 9xl0~ 9 
(CuH), 1.6xl0~ 6 (AlH),4xlO- 8 ( MgH), 3xl0~ 7 (CaH), 2xl0~ 9 
(LiH), 3xl0~ 9 (KH), and 10~ 9 (FeH). We assumed that CaH and 
FeH have a dipole moment of one Debye. In addition, SiH has 
its rotational transitions in the SPIRE and PACS frequency do- 
mains. However, its dipole moment is rather low and the upper 
limit to its abundance relative to H2 is 6xl0~ 6 . Another inter- 
esting species is HF for which our chemical models predict an 
abundance ~10 times below that of HC1. Unfortunately, its J=l-0 
and J=2-l lines are strongly blended with HCN vibrationally ex- 
cited lines and those of other species. The data used to estimate 
these upper limits are shown in Fig. 1 of lDecin et al. I (1201 Oh . 

4. Conclusions 

The detection of HC1 in IRC+10216 indicates the active role 
of chlorine in the chemistry of the warm innermost region of 
the circumstellar envelope. Together with A1C1, HC1 is the most 
abundant chlorine-bearing species in this circumstellar envelope, 
in contrast NaCl and KC1 having abundances 10 and 30 times 
lower, respectively. HIFI observations are needed to spectro- 
scopically resolve the rotational lines of light species and to dis- 
tinguish them from those of more abundant molecules. These 
observations will provide more precise upper limits and might 
perhaps hold some detections of light species. 
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Fig. 1. Continuum-removed spectra of IRC +10216 observed 
with SPIRE and PACS showing the first 6 rotational lines of HCI 
(black histograms) and the line profiles resulting from the LVG 
model (continuous red lines). The contribution of the isotopes of 
CS and the vibrationally excited states of HCN to some of the 
HCI and H 37 C1 lines has been estimated from adjacent lines of 
these species (see text). 




Fig. 2. Abundances of several hydrides computed by thermo- 
chemical equilibrium for the stellar atmosphere and inner en- 
velope of IRC +10216. Abundances are relative to total H and 
are shown as a function of the distance to the center of the star. 



